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Abstract. NLO QCD corrections to Vector Boson Pair Production via Vector 
OO Boson Fusion have recently been calculated and implemented in a parton- 

level Monte-Carlo program with full experimental cuts. We briefly sketch the 
elements of the calculation and show numerical results for the Large Hadron 
Ph Collider. 

^ Introduction 

I— I The Vec;tor Boson Fusion (VBF) process qq qqH is one of the most promis- 

^ ing channels for the discovery of the Higgs particle and the measurement of 

^ its properties at the Large Hadron Collider (LHC) [1]. It proceeds through a 

channel scattering of the two initial-state quarks mediated by a weak bo- 
son, with the Higgs emitted oS^ the boson propagator. The kinematic features 
*0 that make this process phenomenologically relevant are: the presence of two 

highly-energetic jets in the final state, the large rapidity interval between these 
<^ jets and the absence of noticeable jet activity inside this rapidity interval. 

Moreover the Next-to-Leading Order (NLO) QCD corrections to the total 
cross sections [2] and for the differential distributions [3] have been found to 
be quite modest (enhancing the Leading Order (LO) result by 5-10%), thus 
^ ^ pointing towards a good stability of the perturbative result. 

Even though the cross section for the VBF process is somewhat smaller 
;h than the one for the gluon fusion gg ^ H channel at the LHC, the distinctive 

^ features cited above greatly help in distinguishing the signal from the back- 

grounds and make VBF an ideal candidate for Higgs discovery and precision 
measurements. 

One of the most relevant backgrounds to a VBF H — > VV signal is the 
process qq — > qqVV, i.e. vector boson pair production via VBF [4]. It shows 
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exactly the same kinematical features as VBF Higgs production, thus being 
an irreducible background. In addition, it is known that the scattering of lon- 
gitudinal vector bosons is intimately related to the mechanism of electroweak 
symmetry breaking (EWSB), and an enhancement of qq — > qqVV over the 
Standard Model predictions at high center-of-mass energies could be a possible 
signal of strong EWSB (see, for instance, [5] and references therein). 

It is thus clear that an accurate prediction for the electroweak production 
of a vector boson pair plus two jets is mandatory for new physics searches at 
the LHC. 

In the following we will present the results obtained in three recent pa- 
pers where we computed the NLO QCD corrections to the processes qq 
qqW+W- [6], qq qqZZ [7], qq qqW^Z [8], including the full leptonic 
decays of the vector bosons. The calculations have been implemented in a 
fully-flexible parton level Monte-Carlo program, allowing for the computation 
of jet observables and a straightforward implementation of experimental cuts. 

Selected topics from the calculation 

The main challenges of the calculation were the huge number of Feynman 
diagrams involved and the numerical instabilities arising from pentagon con- 
tributions to the virtual part of the cross section. 

When computing a multi-parton process (in this case, a 2—^6 process con- 
sidering full leptonic decays of the vector bosons) one has to find an efiicient 
way to speed up the numerical evaluation. In our Monte-Carlo code the sum 
of sub-amplitudes encountered in several diagrams, and involving only EW 
bosons and leptons, is computed only once per phase-space point (for de- 
tails, see [6]): this method is particularly efficient in the computation of real 
emission corrections, due to the large number of contributing diagrams. 

The soft and coUinear singularities of the NLO real contributions have 
been taken into account by means of the standard Catani-Seymour dipole 
subtraction [9]. Since the divergences only depend on the colour structure of 
the external partons, the subtraction terms are identical in form to the ones 
appearing in Higgs production via VBF [3] . 

The EW bosons exchanged in the channel are colour-singlets, thus there 
cannot be any virtual contribution at 0{as) from gluons attached both to the 
upper and lower quark lines. This allows us to consider virtual radiative cor- 
rections separately for the single quark lines, leading to diagrams containing 
loops with up to five external legs {pentagons). After the cancellation of the in- 
frared poles between real and virtual contributions, the finite remainder of the 
loop amplitudes can be computed by means of the Passariiio-Veltman reduc- 
tion formalism [10] in the case of two-, three- and four-point tensor integrals. 
In the case of pentagons, numerical instabilities show up when kinematical 
invariants, such as the Gram determinant, bec;ome small for some regions of 
phase space. For the pentagons, we have thus used the recently proposed al- 
ternative reduction formalism by Denner and Dittmaier [11], which results in 
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a fraction of numerically unstable events at the per-mille level (for details, see 
[8]). 

Numerical results 

In the following we will present numerical results at NLO QCD accuracy 
obtained with our Monte-Carlo code for EW +W-Wjj,ZZjj,W^Zjj pro- 
duction at the LHC. 

We used the CTEQ6M parton distributions with as^O.llS at NLO and 
the CTEQ6L1 set at LO [12]. We chose toz=9L188 GeV, mw=80A19 GeV 
and Gf=1.166x 10~^ as EW input parameters, obtaining aQ£;£)=l/132.54 
and sin^ 6'vi/=0. 22217. We have set fermion masses to zero, neglecting external 
bottom and top quark contributions. Jets have been reconstructed by means 
of the fc^-algorithm [13, 14] with resolution parameter D—0.8. 

Typical VBF cuts have been imposed: here we show those used in the 
VV^Z case: 

• two hard "tagging" jets: ptj >20 GeV, 
|yj|<4.5, Mjj > 600 GeV 

• large rapidity separation between jets: '^yjj >4, yji x yj2 <0 

• leptoncuts: pr; >20 GeV, |r?;| <2.5, m,/ >15 GeV 

• "separation" cuts: ARji >0.4, ARa >0.2 , 

where ARji and ARu denote the jet-lepton and lepton-lepton separation 
in the rapidity-azimuthal angle plane, respectively, and mu the invariant mass 
of an electron or muon pair. 
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Fig. 1. Scale dependence of the total cross section for W^Z production via VBF 
for two values of the central scale (see text). 
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In Figure 1 (from [8] ) we show the total cross section displayed as a function 
of the factorization and renormalization scales fJ.F,R = ^F.RfJ-o the case of 
W^Zjj production. We have considered two possible values for the central 
scale: jj,o = {mz + mw)/'^ or = Q, where Q is the momentum transfer 
carried by the exchanged vector boson in VBF graphs (for details, see [8]). 
The K-factor is K=Q.%1 in the first case and iir=1.04 in the second one. In 
both cases the scale dependence, for instance in the range 0.5 < ^ < 2, is 
greatly reduced when passing from LO (about 10%, dotted black curves) to 
NLO (about 2%). At NLO, we have considered three cases: = = £, (solid 
red lines), ^f — ^,^r — ^ (dashed green lines), ^f = ^,^r = 1 (dot-dashed 
blue lines). 
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Fig. 2. Transverse-momentum distribution of the tagging jet with the highest px 
in W^W' production via VBF (left), with the corresponding if- factor (right). The 
cuts of Ref. [6] are used. 



In Figure 2 (from [6]) we present the transverse-momentum distribution 
of the tagging jet with the highest px in the WWjj case, together with the 
corresponding K-factor. The figure shows a strong change in shape when going 
from LO to NLO, with a 10-20% enhancement of the cross section at low values 
of pt {pt <100 GeV) and a corresponding decrease at higher px values: this 
effect is mainly due to the extra parton coming from real emission at NLO. 

Finally, in Figure 3 (from [7]) we show the differential distribution for EW 
ZZjj production with respect to the invariant mass AIzz without (left) and 
with (right) the inclusion of the Higgs contribution for a scale fi — Q. Apart 
from the pronounced resonance behaviour visible in the right plot, we note 
that the LO and NLO predictions are virtually indistinguishable in both cases, 
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Fig. 3. Distribution of the invariant mass Mzz in ZZ production via VBF without 
(left) and with (right) the Higgs boson contribution. The cuts of Ref. [7] are used. 



indicating an excellent stability of the perturbative calculation for this scale 
choice. 

Acknowledgements. The work of B.J. is supported by the Japan Society 
for the Promotion of Science. The work of G.B. is supported by the Deutsche 
Forschungsgemeinschaft under SFB TR-9 "Computergestiitzte Theoretische 
Teilchenphysik" . 

References 

1. D. Zeppenfeld, R. Kinnunen, A. Nikitenko, and E. Richter-Was, Phys. Rev. 
D 62, 013009 (2000) [arXiv:hep-ph/0002036]; D. Zeppenfeld, in Proc. of the 
APS/DPF/DPB Summer Study on the Future of Particle Phystcs, Snowmass, 
2001 edited by N. Graf, eConf C010630, p. 123 (2001) [arXiv:hep-ph/0203123]; 
A. Belyaev and L. Reina, JHEP 0208, 041 (2002) [arXiv:hep-ph/0205270]; 
M. Diihrssen et al., Phys. Rev. D 70, 113009 (2004) [arXiv:hep-ph/0406323]. 

2. T. Han, G. Valencia and S. Willenbrock, Phys. Rev. Lett. 69, 3274 (1992) 
[arXiv:hep-ph/9206246] . 

3. T. Figy, C. Oleari and D. Zeppenfeld, Phys. Rev. D 68, 073005 (2003) 
[arXiv:hep-ph/0306109] . 

4. D. Rainwater and D. Zeppenfeld, Phys. Rev. D 60, 113004 (1999) [Erratum- 
ibid. D 61, 099901 (2000)] [arXiv:hep-ph/9906218]; N. Kauer, T. Plehn, 
D. Rainwater, and D. Zeppenfeld, Phys. Lett. B 503, 113 (2001) [arXiv:hep- 
ph/0012351]. 

5. M. S. Chanowitz, Czech. J. Phys. 55, B45 (2005) [arXiv:hep-ph/0412203]. 

6. B. Jager, C. Oleari and D. Zeppenfeld, JHEP 0607, 015 (2006) [arXiv:hep- 
ph/0603177]. 



6 



Giuseppe Bozzi, Barbara Jager, Carlo Oleari and Dieter Zeppenfeld 



7. B. Jager, C. Oleari and D. Zeppenfeld, Phys. Rev. D 73, 113006 (2006) 
[arXiv:hcp-ph/0604200]. 

8. G. Bozzi, B. Jager, C. Oleari and D. Zeppenfeld, Phys. Rev. D 75, 073004 
(2007) [arXiv:hep-ph/0701105]. 

9. S. Catani and M. H. Seymour, Nucl. Phys. B 485, 291 (1997) [Erratum-ibid. 
B 510, 503 (1998)] [arXiv:hep-ph/9605323]. 

10. G. Passarino and M. J. G. Vcltman, Nucl. Phys. B 160, 151 (1979). 

11. A. Denner and S. Dittmaier, Nucl. Phys. B 658, 175 (2003) [arXiv:hep- 
ph/0212259], Nucl. Phys. B 734, 62 (2006) [arXiv:hep-ph/0509141]. 

12. J. Pumplin, D. R. Stump, J. Huston, H. L. Lai, P. Nadolsky and W. K. Tung, 
JHEP 0207, 012 (2002) [arXiv:hep-ph/0201195]. 

13. S. Catani, Yu. L. Dokshitzcr, and B. R. Webber, Phys. Lett. B 285 291 (1992); 
S. Catani, Yu. L. Dokshitzer, M. H. Seymour, and B. R. Webber, Nucl. Phys. 
B406 187 (1993); S. D. Ellis and D. E. Soper, Phys. Rev. D 48 3160 (1993). 

14. G. C. Blazey et al., arXiv:hep-ex/0005012. 



